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Abstract 

The restricted planar elliptic three body problem models the motion of a massless body 
under the Newtonian gravitational force of two other bodies, the primaries, which evolve in 
Keplerian ellipses. 

A trajectory is called oscillatory if it leaves every bounded region but returns infinitely 
often to some fixed bounded region. We prove the existence of such type of trajectories for 
any values for the masses of the primaries provided the eccentricity of the Keplerian ellipses 
is small. 
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1 Introduction 

The restricted planar elliptic three body problem (RPE3BP from now on) models the motion of 
a body of zero mass under the Newtonian gravitational force of two other bodies, the primaries, 
which evolve in Keplerian ellipses with eccentricity eo- Without loss of generality one can assume 
that the masses of the primaries are p and 1—/i, their period is 27r and their positions are —/-iqo(t) 
and (1 — n)qo(t), where 

q 0 (t) = ( p(t ) cos v(t), p(t) sin v(t)) 

with 


1 + eo cos v(t) 

and v(t) is called the true anomaly, which satisfies u(0) = 0 and 

dv (1 + eo cos v ) 2 
dt = (1 — eg) 3 / 2 ' 

Then, the motion of the third body is described by the following Hamiltonian 

H(q,p,t;e o) = - V{q,t\e 0 ) (1) 

where 

V{q,t,e 0 ) || g + fj, qo {t)\\ + ||g - (1 - fJ,)q Q (t )|| 

where q,p G M 2 . 

The eccentricity of the ellipses satisfies eo G [0,1). If eo = 0, the primaries describe circular 
orbits. This case is known as the restricted planar circular three body problem (RPC3BP from 
now on). In this paper we consider eo > 0 and small. Note that the Hamiltonian also depends 
on fi G [0,1/2]. We do not write this dependence explicitly since for us p is a fixed positive 
parameter. 

The purpose of this paper is to analyze some particular orbits of this Hamiltonian system: 
the oscillatory motions, that is, orbits which leave every bounded region but return infinitely 
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often to some fixed bounded region. To prove the existence of such orbits, we use the framework 
usually considered in the study of Arnold diffusion in nearly integrable Hamiltonian systems. 

If one takes fi = 0, the system reduces to a central force problem which is integrable and 
therefore cannot have oscillatory motions. That is, they can take place for the three body 
problem but not in its limit // —»• 0. In this paper we prove that oscillatory motions are possible 
for any value of /i £ (0,1/2] and eo > 0 small enough. The existence of such motions for any 
(i and eo = 0 (the circular problem) has been recently proved by the authors of this paper in 
|CMSl5| . 

To obtain oscillatory orbits we work in a nearly integrable setting and we use perturbative 
methods that provide invariant objects which can be analyzed and act as a skeleton that such 
orbits follow. Nevertheless, since we give results for any /r € (0,1/2], the nearly integrable setting 
cannot be in terms of /r —> 0. Instead, we consider the following regime. Take the body of zero 
mass very far away from the two primaries. Then, at first order the third body perceives the two 
primaries as just one body at the center of mass plus a periodic perturbation whose smallness 
comes from the ratio of the distance between the two primaries over the distance between the 
third body and the center of mass. Then, taking the ratio of distances small enough, one has 
an integrable Hamiltonian plus a small periodic perturbation (see jGMS15j ). 

1.1 Final motions in the three body problem 

One of the most important questions in the analysis of the three body problem (either restricted 
or non restricted, planar or spatial) is the study of the final motions. That is, what type of 
behaviors can happen as time t —> Too. Its analysis was initiated by Chazy in 1922, when he 
gave a complete classification of the possible final motions (see Section 2.4 of [ AKN88I 1. In the 
restricted setting, the possible final motions are the following: 

• H ± (hyperbolic): \\q(t)\\ -T oo and ||(j(t)|| —> c > 0 as t —> Too. 

• P ± (parabolic): \\q(t)\\ -T oo and ||(j(t)|| —> 0 as t —> Too. 

• B ^ (bounded): limsup t ^ ±00 ||g|| < Too. 

• OS^ (oscillatory): limsup^^^ ||g|| = Too and liminf^j-oo ||(?|| < Too. 

Examples of all these behaviors, except the oscillatory motions, were already known by Chazy. 
The other three behaviors certainly already exist for the two body problem, where motion 
is confined to conics in the state space and hyperbolic motions arises on hyperbolas, parabolic 
motions on parabolas and bounded motion on ellipses. As already explained, oscillatory motions 
cannot appear in the two body problem. 

The study of oscillatory motions was initiated by Sitnikov in the sixties jSitfiOj . He proved 
their existence in a very carefully chosen symmetrical model consisting of two bodies of equal 
mass revolving in planar ellipses around their center of mass and a third body of mass zero 
moving along the perpendicular axis at the center of mass. His work also proved that X~f]Y + ^ 
0 for X,Y = H, P, B, OS. 

Later, Moser [Mos73j gave a new proof of the Sitnikov result. The subsequent results in 
the area strongly rely on the ideas developed by Moser. The present paper certainly also relies 
on some of them but also uses other ideas recently developed for the study of Arnold diffusion. 
Using Moser ideas, Llibre and Simo |LS80b| obtained oscillatory motions for the collinear three 
body problem. 
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For the planar three body problem, the first result was also by Llibre and Simo |LS80a| , where 
they proved the existence of oscillatory motions for the RPC3BP for small enough values of the 
mass ratio p (Hamiltonian dTJ) with eo = 0). Using a real-analyticity argument, Xia |Xia92j 
extended their result to any p G (0,1/2] except a finite (unknown) number of values (Moser 
|Mos73| had previously used this argument in the Sitnikov problem). The first general result 
was obtained by three of the authors of the present paper [GMS15] , who proved the existence 
of oscillatory motions of the restricted planar circular three body problem for any value of the 
mass ratio /j, G (0,1/2]. 

A completely different approach using Aubry-Mather theory and semi-infinite regions of 
instability was developed in [GKlll IGKlObl IGKIOaj . The authors considered the RPC3BP 
with a realistic mass ratio for the Sun-Jupiter system. Using computer-assisted methods, they 
proved the existence of orbits with initial conditions in the range of our Solar System which 
become oscillatory as time tends to infinity. 

All the mentioned results deal with systems which can be reduced to a two dimensional 
area preserving map. This makes the proof of existence of oscillatory motions considerably 
simpler than if one considers other three body problems (restricted elliptic, restricted spatial, 
non-restricted) which have more degrees of freedom. For higher dimensional systems, Moser 
ideas are harder to apply and the results are more scarce. 

In more degrees of freedom the first result is due to Alexeev |Ale69| (actually before Moser 
ideas) who generalized Sitnikov result [Sit60| to small mass for the third body. Positive mass 
increases the number of degrees of freedom of the system to three. Moeckel |Moe84l IMoe07| has 
proven results concerning oscillatory motions in the regime of small angular momentum (and 
thus close to triple collision). 

Concerning the RPE3BP (pQ), Robinson gave a conditional result on the existence of oscilla¬ 
tory orbits jRob841 lRob!5j . He proved their existence provided certain homoclinic points were 
present. As far as the authors know, the existence of such points has not been proved yet. 

A fundamental problem concerning oscillatory motions is to measure how abundant they are. 
As pointed out in CK 12 , in the conference in honor of the 70th anniversary of Alexeev, Arnold 
posed the following question: Is the Lebesgue measure of the set of oscillatory motions positive? 
Arnold considered this problem the central problem of celestial mechanics. Alexeev conjectured 
in |Ale71j that the Lebesgue measure is zero (in the English version |Ale81j he attributes this 
conjecture to Kolmogorov). This conjecture remains wide open. 

The only result dealing with the abundance of oscillatory motions is the recent paper |GK12| . 
The authors study the Hausdorff dimension of the set of oscillatory motions for the Sitnikov 
example and the RPC3BP. Using |Mos73] and |LS80aj . they prove that the Hausdorff dimension 
of the set of oscillatory rnotins is maximal for a Baire generic subset of an open set of parameters 
(the eccentricity of the primaries in the Sitnikov example and the mass ratio and the Jacobi 
constant in the RPC3BP). 

1.2 Arnold diffusion in the three body problem and growth in angular mo¬ 
mentum 

As we have said, we study a regime where the RPE3BP is close to a two body problem. In the 
latter, the angular momentum is a first integral. This fact is no longer true in the former. Thus, 
a natural question is whether the angular momentum of the zero mass body only varies by a 
small amount (with respect to the perturbative parameter) or can make big excursions. 

This question fits into the framework of Arnold diffusion }Arn64] . Consider a nearly inte- 
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grable system in action-angle coordinates 

H(tp,I) = ifoCO + £T#i (¥>,/), <P £ r, I £V C l n , e « 1. 

For e = 0, the action variables / are constants of motion. Arnold diffusion analyzes the drastic 
changes that the actions can undergo for small e > 0 . 

In the setting of nearly integrable IV-body problems, the existence of Arnold diffusion can 
be analyzed in very different regions in the phase space and for very different ranges of the 
involved parameters. As far as the authors know, the first paper dealing with Arnold diffusion 
in Celestial Mechanics is |Moe96| (later completed in [ZhelOj ) who considers the five body 
problem. In [DGRllj . the authors analyze unstable behavior for the three body problem close 
to the Lagrangian point L\. 

Concerning to the growth of angular momentum, the paper [FGKRll] obtains such behavior 
for the RPE3BP (Hamiltonian J!Q)) along the mean motion resonances, which implies a change 
of eccentricty in the osculating ellipse of the body of mass zero. More related to our paper, is 
the recent |DKdlR,S14j . where such behavior is obtained in a neighborhood of parabolic motions, 
the so-called invariant manifolds of infinity. In |DKdlRS14| . it is proven the existence of orbits 
whose angular momentum G(t) satisfies G(0) < G i, G(T) > G 2 for some T > 0 and for any 
given G *2 > G\. Nevertheless, they need to assume that p is exponentially small with respect to 
G 2 and eo is polynomially small with respect to G 2 . 

We expect that this result can be generalized to any value of the mass ratio. Nevertheless, 
in the present paper we only give a conditional result (see Remark m- 

1.3 Main result 

We give in this section the main result that we obtain. 

Theorem 1.1. Fix any p £ (0,1/2]. There exists e${p) > 0 such that for any eo £ (0,ej^//)) 
there exists an orbit ( q(t),p(t )) of (JTJ) which is oscillatory. Namely, it satisfies 

limsup \\q\\ = +00 and liminf ||g|| < + 00 . 

t— H -00 t-H-oo 

This theorem shows that OS + 7 ^ 0. Proceeding analogously, one can show that OS~ 7 ^ 0. 
Nevertheless, our techniques do not allow us to show that OS~ H OS + 7 ^ 0. To do that, one 
needs to consider more sophisticated techniques. Indeed, the ideas developed by Moser rely on 
the construction of a horseshoe which has branches arbitrarily close to infinity. This allows to 
combine all possible past and future behavior. Here, as we explain in Section 11.41 we rely on 
different techniques, which are simpler to generalize to higher dimensions: shadowing certain 
invariant objects. These techniques can only be applied to one time direction, either future or 
past, but not to both of them at the same time. 

We only obtain oscillatory motions provided the primaries perform nearly circular orbits. It 
is expected that such motions exist for any value of the eccentricity eo £ [0,1). The approach 
presented in this paper can be applied to this more general setting. The only additional difficulty 
is to extend the proof of transversality of the invariant manifolds of infinity to this wider range 
of parameters (see Section [2] for more details). In this paper, we use perturbative arguments 
and use that this transversality is already known for the circular problem }GMS15j . 

Remark 1.2. The geometric framework we use to prove Theorem \l.l\ can be also used to obtain 
orbits with large drift in the angular momentum, provided some non-degeneracy condition is 
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satisfied. If such condition could be verified to be true, one would generalize the results in 
IDKdlRSlIfl by obtaining the following. Fix fi € (0,1/2]. There exists Go(/T) > 0 such that for 
any G 2 > G\ > Go(n), there exists eQ(G 2 ) > 0 such that for any eo € (0, eo(G 2 )) there exist a 
time T > 0 and an orbit ( q(t),p(t )) of (HD whose angular momentum G satisfies 

G(0) < G\ and G(T) > G 2 . 

The needed non-degeneracy condition is explained in R,emark \2.6l 

1.4 Common framework for oscillatory motions and Arnold diffusion 

The purpose of this section is to relate the results of both Theorem 11.11 and Remark 11.21 and put 
them in the same framework. This is explained in more detail in Section [2] 

Let us start by explaining Moser ideas to obtain oscillatory orbits for the Sitnikov example 
[Mos73] and applied in [LS80a ] to the RPC3BP. The RPC3BP (Hamiltonian (pQ) with eo = 0) 
has a first integral, the Jacobi constant. Then, in suitable coordinates, the RPC3BP can be 
reduced to a two dimensional Poincare map for which “infinity” {|g| = + 00 , q = 0} is a parabolic 
fixed point. 

Assume for a moment that this fixed point is hyperbolic. Then, it would have stable and un¬ 
stable invariant manifolds. Assume that these invariant manifolds intersect transversally. Then, 
Smale Theorem would ensure that there exists a horseshoe. Suitable orbits in this horseshoe 
travel close to the invariant manifolds and the lim inf of the distances to the fixed point is zero. 
Such orbits, in the original coordinates, are oscillatory. 

The infinity point of the RPC3BP is not hyperbolic but parabolic. |LS8flaj . as done in 
(Mos73| for the Sitnikov problem, shows that even if it is parabolic one can carry out the same 
strategy. Note that one has to face different difficult issues: prove the existence of the invariant 
manifolds of the parabolic point (see |McG73] ), prove that they intersect transversally and prove 
a Lambda lemma for parabolic points which implies the existence of symbolic dynamics. 

To carry out this strategy in the elliptic case is certainly more involved since the phase space 
has dimension five and therefore one cannot reduce the dynamics to an area preserving map. 
The stroboscopic Poincare map is four dimensional. Thus, infinity cannot be reduced to a fixed 
point but it forms a cylinder with one angular variable and one real variable (see Section ED- 
This cylinder is “normally parabolic” and it has invariant manifolds. Since in this paper we 
consider the elliptic problem as a perturbation of the circular one, the result in [GMS15 ] implies 
that the invariant manifolds of this cylinder intersect transversally. 

To prove that such transversal intersections lead to oscillatory motions we do not rely on the 
construction of symbolic dynamics as Moser did (Mos73j . Instead, we use the standard method 
of Arnold diffusion to construct a transition chain of tori [Arn64| (in the present setting of fixed 
points). That is, we find a sequence of fixed points belonging to the cylinder of infinity which are 
connected by transversal heteroclinic orbits. In Arnold diffusion problems usually it is enough 
to construct a finite chain. Instead, to obtain oscillatory orbits one has to construct an infinite 
chain and then prove the existence of an orbit which shadows the chain. This construction is 
much simpler than a horseshoe and leads to oscillatory motions forward (or backward) in time. 
Certainly a horseshoe gives much more information and imply a plethora of different types of 
motion. In particular, it allows to combine different types of final motions in the past and in the 
future (including motions which are oscillatory as both time tends to plus and minus infinity). 
Our simpler techniques can only show the existence of oscillatory motions in one direction of 
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time. Nevertheless, we think that one of the main interests of this paper is to show that such a 
simple mechanism as the existence of transtion chains leads to oscillatory motions. 

The vertical direction of the cylinder of infinity can be parameterized by the angular mo¬ 
mentum of the orbits. To prove the existence of oscillatory orbits one needs to construct a 
chain which remains in a compact portion of this cylinder. The reason is that orbits shadowing 
chains with unbounded angular momentum cannot be oscillatory since liminf|g| = +oo. In 
fact, we can prove the existence of chains confined in very thin portions of the cylinder and 
thus with almost constant angular momentum. This is done by considering the scattering map 
I fill I ..SO* and studing its dynamical properties. Now, one can ask the opposite question. It is 
possible to construct a chain which implies a large deviation in the angular momentum? Orbits 
shadowing such chains have a big drift in angular momentum and present the phenomenon of 
Arnold diffusion as explained in Remark 11.21 Summarizing, the orbits given both by Theorem 
o and Remark 11,21 are obtained thanks to suitable transition chains associated to the cylinder 
of infinity. 

The structure of the paper goes as follows. First in Section [2] we analyze the invariant 
manifolds of infinity. We prove its existence and regularity and we analyze their transversal 
intersections. This allows us to construct the needed transition chain of periodic orbits. In 
Section [3l we state a Lambda lemma which can be applied to the invariant manifolds of the 
normally parabolic cylinder of infinity. This Lambda lemma allows us to prove Theorem 11.11 
The prove of the Lambda Lemma is deferred to Section 01 


2 The invariant manifolds of infinity 


Let (r, a) be the polar coordinates in the plane and (y, G) their symplectic conjugate momenta. 
Then, the Hamiltonian ([I]) becomes 


H(r,a, y, G, f; e 0 ) = 


1 / G 2 


+ y‘ 


U(r,a,t;e 0 ), 


where U ( r , a , t; eo) = V ( re ia , t; eo). The variable G is the angular momentum of the third body. 
As we deal with a non-autonomous system, we add the equation s = 1 to the Hamiltonian 
equations of Hamiltonian H(r, a, y, G, s; eo). 

Since we want to study the invariant manifolds of infinity, we consider the McGehee coordi¬ 
nates ( x,a,y,G) where 

r = — 7 T, for x > 0. 
x z 

We obtain the new system 


x = 

a = 

s = 


;x y 


y = -G 2 x 6 
y 8 

dU 


G G = 8a 


x 3 8U 

4 dx 


where the potential U is given by 


U(x,a,s\e o) = C/(2/x , a, s;eo) = — 


x 2 (1 — y, n 


H- 

^5 


( 2 ) 
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with 


\q - <?s|" = 0 's = 1 - h r o {s)x z cos (a - v(s)) + -/Tro(s)a ; 4 


\q-qj\ = Cj = 1 + (1 - n)ro{s)x cos (a - v(s)) + -(!-//) r 0 (s)x . 


We write the potential as 


Since 


U = — + AG. 

U(r,a,t;e 0 ) = - - - 3cos(a - v(t))Z-^- + O (ZZ\ , 

r z r° \r*/ 


( 3 ) 


r° \r* 

the potential A U satisfies A U = 0(px 6 ). 

In view of (|2|). now “infinity” is foliated by the parabolic periodic orbits 

^■a 0 ,G 0 = {(x, a, y, G, s)gMxTxMxRxI|x = |/ = 0, a = ao, G = Go}. 


Next theorem claims that each periodic orbit A a 0 j G 0 has stable and unstable 2-dinrensional 
invariant manifolds W s (A aQ! G 0 ) and W u {h. ao ^G 0 ). 

Theorem 2.1. Let fa the flow of the system @ and define the projections n x (x. a, y, G,s) = x 
and Tr^ Xjy )(x,a,y,G,s) = (x,y). Let (ao,Go) € T x R. There exists po such that for any 
0 < p < po, the local stable set 


W s p (A ao ,Go) = {( x,a,y,G,s ) eKxTxKxKxI) Tr x fa(x, a, y,G, s) > 0, 

W(x, y )fa{x,a,y,G,s)\ < p, t Kmdist(fa(x,a,y,G,s),A aO:Go ) = 0}, 
is a 2-dimensional manifold. 

Moreover, there exists uq > 0 such that Wp(A ao> G 0 ) admits a C°° parametrization Gq : 
[0, rto) xT->KxIxlxlxT, with 7 a 0 g 0 (0, s) = (0, ao, 0, Go, s), analytic in (0,«o) x T, 
which depends analytically on (ao, Go) G T x M. The analogous result for the (local) unstable set 
also holds. 

The proof of this theorem follows from Proposition 14.11 in Section [4] (the technical details 
are done in Appendix O. A related result concerning the invariant manifolds of infinity is given 
in [Rob84j (see also |Robl5] h Nevertheless we have included the proof of this theorem, which 
is based in the parameterization method [CFdlL03l !BFdlLM07j . for two reasons. One the one 
hand, we are not able to follow all the details in the proof in jE.ob84llRobl5| . On the other hand, 
our result provides the regularity of the stable/unstable foliations of the invariant manifolds, 
needed in Lemma 13.11 

Theorem 12.II shows that the points which tend asymptotically in forward (or backward) time 
to the periodic orbit Aq, 0 i g 0 form a manifold. Since this periodic orbit is not hyperbolic but 
parabolic the rate of convergence is not exponential and its invariant manifolds are not analytic 
at the periodic orbit but only C°°. They are analytic at any other point. 

As we have explained in Section [I] when p = 0, the invariant manifolds l / F s (A Q 0 ] c' 0 ) an d 
W u (A a0tGo ) coincide and form a two-parameter family of parabolas in the configuration space 
(in the original cartesian coordinates). We study the splitting of these invariant manifolds for 
any p G (0,1/2] in two steps. First for the circular problem (eo = 0) and then for the elliptic 
problem (0 < eo <S 1). In this study we need formulas for the homoclinic manifolds when p = 0. 
For the derivation of such formulas, one can see |LS 8 flaj . 
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Lemma 2.2. Take fj, = 0 and fix £ T and Go 0. Then, system ((21) has a family of 
homoclinic orbits to the periodic orbit A QOt c 0 , which is given by 

Xh ^ Go) = Go (l + r 2 ) l/2 

yh ^ ;Go) = G 0 (l + r 2 ) 

ah(^;ao) =ao + 5h (t), 5h(i) = 2arctanr 

G h (t;G 0 ) =G 0 

s h(^j s o) =s 0 + t 

where sq £ T is a free parameter and r and the time t are related through 



From now on, we abuse notation and we consider the functions Xh, y h and a h both as 
functions of r or t. 


2.1 Scattering map for the circular problem in the extended phase space 

Using the results of Theorem 12.11 for any G\ > 0, the invariant set 

A [Gl ’ +oo) = \J A QOiGo = {(x, a, y, G, s) £ R x T x R x R x T | x = y = 0, G > Gi} (4) 
«oGT,Go>Gi 

is a “normally parabolic” 3-dimensional invariant manifold with stable and unstable 4-dimensional 
invariant manifolds 


^(AlGi’+oo)) = U 1 U ? (A Qo ,Go), S = u,s. 

ao£T,Go>Gi 

Theorem 2.2 of [ QMS 15] implies that, when eo = 0 and p, £ (0,1/2], there exists G* 1 
such that, for any G\ > G* , the invariant manifolds W s (A[ Gl,+0 °)) and W u (A[ G i,+00 )) intersect 
transversally in the whole space along two different 3-dimensional homoclinic manifolds. 

Following [DdlLS08l . this transversality allows us to define two scattering maps 5^ as¬ 
sociated to the two different transversal homoclinic intersections between W s (A[ Gl,+0 °)) and 
W“(A[ Gi,+0 °)) and to obtain formulas for these maps. To this end, we consider the Poincare 
function for eo = 0, 

/ CO 

AU(x h (a + t; G 0 ), a 0 + a h (a + t; G 0 ), s 0 + t\ 0) dt (5) 

-OO 

where A U is the potential defined in ((3]) and (xh, 5h) are components of the parameterization 
of the unperturbed separatrix given in Lemma 12.21 

For any (cto, Go, so)> the function a i->- L(ao, Go, so, er; 0) has two critical points erf given by: 

erf = so - ao, <x* = 7r + s 0 ~ Go- (6) 

This fact is given by Proposition 3.1 of (GMS15| (note that this proposition is stated in certain 
scaled variables). 


9 










The reason to obtain such simple formulas for the critical points is twofold. On the one 
hand, when eo = 0, the potential only depends on the angles ao and so through ao ~ so ■ Thus, 
since 

L(ao, Go, s 0 , cr; 0) = L(a o, Go, so — u, 0; 0), 

the Poincare function only depends on one angular variable ao — so + a. On the other hand, its 
Fourier expansion only contains cosines (see the Appendix of | iGMS15j l. 

Associated to the zero <j 5_, there is a heteroclinic connection between two periodic orbits in 
A[ g 1’ + °o), which are /rG^-close to A Q0 ! g 0 (analogously for a *_). These heteroclinic connections 
satisfy 

r aoTo,s 0 ( t ) _ Go), 2/h(o± + t ; Go), ao + ah(o± + 1 ; Go), Go, so + 1) + O (/iG 0 4 ) 

(see [GMS15] ). Note that these heteroclinic orbits are well defined for any ao € T, so G T and 
Go > Gi (see (|T|)). This implies that there is a homoclinic channel (see jDdlLSOBj l which is 
defined for all points in Ai Gl,+00 ). Thus, we can define global scattering maps 

£± . ^\Gi,+oo) _^ ^[Gi,+oo) 

Following [ DdlLSOB] . recall that x + = 5^ (x-) if there exists a heteroclinic connection between 
these two points through the prescribed homoclinic channel (see |DdlLS08] for a more precise 
definition and properties of the scattering map). Usually it is not possible to define globally the 
scattering map since it is only defined locally in open sets. Then, globally, it can be multivaluated 
(see [DdlLSOfi] !. The particular form of the circular problem allows us to define it globally. Note 
also that, in principle, the scattering map should map A^ 1 ’ 00 ) to a bigger cylinder. Nevertheless, 
as it is shown in the next proposition, in the circular case the image is the same cylinder. 

Proposition 2.3. LetG\>G*. The scattering maps S^ 1 : Tx [Gi,+oo) xT — »Tx[Gi,+oo)x 
T are of the following form, 

< S±(a,G, S ) = (a + / ± (G),G,s), 

where 

q-t r 

/ ± (G) = -Ml-/i)^4+0(G- 8 ). (7) 

This proposition is proven in Appendix [Bj 

2.2 Reduction to the Poincare map 

We reduce the dimension of the system by considering the stroboscopic-Poincare map associated 
to the section S = {s = sq}, 


V : £ 

(x,y,a,G) 



£ 

' P(x,y,a,G ) 


( 8 ) 


Then A Q0 i g 0 = A QQ y; 0 n £ is a two parameter family of parabolic fixed points of V with 1- 
dimensional stable and unstable manifolds 


^(A QOiGo ) = w q { A« 0 ,g 0 ) n £, g = u,s. 
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Analogously, A[ Gl ’+°°) = Al Gl,+00 )nS is the 2-dimensional normally parabolic invariant cylinder 
of infinity with 3-dimensional invariant stable and unstable manifolds 

( A [Gi,+°o)) = ^ a [Gi >+ oo)) n Sj ? = u> 

which, for eo = 0, intersect transversally along two 2-dimensional homoclinic channels. The two 
scattering maps associated to these homoclinic channels are given by 

S±(a,G) = (a + f ± (G),G), (9) 

where is the function given in (J7|) . They do not depend on the section S. 

2.3 Scattering map of the elliptic problem 

Once we have analyzed the splitting of the invariant manifolds of A[ Gl ’ +00 ) for the circular 
problem and derived formulas for the two scattering maps, now we consider the elliptic problem 
for eo small enough. Note that eo is a regular parameter of the elliptic problem and, therefore, we 
can apply classical perturbative arguments to the stroboscopic-Poincare map V (see [ DdlLS08] ). 

Let G *2 > G\ > G* (see ([3]) and Proposition ^. 31) be fixed. We call A^ 1,02 ! to A[ Gl,+00 )n{Gi < 
G < G 2 }, which is compact and invariant. Then, for eo small enough, the stable and unstable 
manifolds of A^ 1 ’^ intersect transversally. Note that the smallness of eo depends on the chosen 
interval. The perturbative arguments imply that there are two global homoclinic channels as in 
the circular problem. These two channels define two scattering maps 

£± . A [Gi,G 2 ] _ s, a [G*,+oo) 


which depend regularly on eo, 


+ eoSf + O (eg) , (10) 

where 5^ are the scattering maps of the circular problem given by Q. The theory developed 
in [DdlLS08] does not directly apply to this case since the invariant cylinder is not hyperbolic, 
but parabolic. Nevertheless, the arguments in |DKdlRS141 Proposition 4] show that the theory 
of scattering maps do apply also to this problem. Hence, the maps are area preserving maps 
on the cylinder. 

To construct oscillatory orbits of the elliptic problem, we need an infinite transition chain: a 
sequence of fixed points belonging to Al Gl with transversal heteroclinic connections between 
consecutive points of the chain. By the definition of the scattering map, any bounded (forward 
or backward) orbit of one of the scattering maps provides such a chain. 

To obtain bounded orbits of the scattering maps we observe that, if eo is small enough, the 
maps S ± possess invariant curves. Indeed, in view of (|T|). are twist maps if eo is small 
enough. Then, we can apply the following twist theorem, due to Herman, from [Her83j : 

Theorem 2.4 (Twist theorem). Let f : [0,1] x T —> [0,1] x T be an exact symplectic C l map 
with l > 4. Assume that f = fo + Sfi, where fo(I,f>) = (L, ip + A(I)), A is C l , \diA\ > M and 
Why < I- Then, ifS^M- 1 = p is sufficiently small, for a set of u of Diophantine numbers 
of exponent 9 = 5/4, we can find invariant curves which are the graph of C l ~ 3 functions u^, the 
motion on them is C l ~ 3 conjugate to the rotation by w, and ll'it^llc'- 3 — C5 1 / 2 . 
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In our setting we have that the twist condition satisfies \dcA\ > gG~ 5 whereas the <5 < eo- 
Therefore, for any G 2 > G\ > G* and eo small enough, the map tS 1 * 1 has KAM curves inside 
A[Gi,G 2 ] Any orbit of iS 1 * 1 in any of these KAM curves is bounded, and so are all orbits between 
any two of these KAM curves. 

Corollary 2.5. Consider any constants G 2 > G\ > G*. Then, for eo small enough (which may 
depend on G 2 and G\), the scattering maps considered in (Hoi) have orbits which remain for 
all time in A( Gl,G2 (. 


As we have explained, coming back to the original Poincare map (j 15 H . this bounded orbit of 
the scattering map corresponds to a sequence of fixed points {Ac^G'j.lfcgN C A( Gl ' G2 ( such that 
W u (h. ak) G k ) intersects transversally iy s (A( Gl,G 2 l) at a point P which belongs to W s {K ak+1) c k+1 )- 


Remark 2.6. To obtain oscillatory motions we have looked for an infinite transition chain with 
bounded angular momentum, or equivalently, for a bounded orbit of one of the scattering maps 
5^. If we want orbits with a drift in the angular momentum (see Remark ]l.'A) . we have to look 
for a transition chain connecting fixed points of V in (USD with a large difference in angular 
momentum. In this case it suffices to look for a finite transition chain. 

For the oscillatory motions, we have used the fact that the scattering maps are nearly 
integrable and possess KAM tori. KAM tori act as barriers for the orbits of the scattering 
maps and therefore, using only one scattering map, it is impossible to construct orbits with a 
large drift in G. The usual strategy to prove the existence of Arnold diffusion (see for instance 
\Arn6 f , \PdlLSOflj ) is to combine one scattering map with the inner dynamics of the invariant 
cylinder induced by the Poincare map (USD. Combining these two dynamics one obtains Arnold 

diffusion. Here this approach is not possible since the inner dynamics is trivial (the cylinder 
A [Gi,G 2 ] 

is filled by fixed points). Thus, we rely on an idea developed in JPKdlRSlIf) which is to 
combine the two scattering maps S + and S~. 

In fLC07\l , Le Calvez showed the following: consider two area preserving twist maps on a 
cylinder. Assume that they do not have common invariant curves. Then, combining them one 
can obtain orbits with a drift in the action component. Thus, as long as the maps S± did not have 
common invariant curves, we could obtain transition chains with a large drift in G. This fact 
has proven to be true for p small enough in fPKdlRSl^ . Nevertheless, it is not straightforward 
to verify it in the present setting. That is the non-degeneracy assumption mentioned in Remark 


3 Shadowing orbits 

The second ingredient in the proof of Theorem II.II is what is usually called a Lambda or Shad¬ 
owing lemma. That is, to analyze how orbits close to the stable manifold of one of the fixed 
points of Al Gl ’ G2 l (see Theorem 12.11) evolve and stretch along the unstable invariant manifold 
of the fixed point. Such results allow to shadow a concatenation of heteroclinic orbits which 
connect different (or the same) fixed points. 

Usually in the literature of Arnol’d diffusion, one considers more general Lambda lemmas 
which allow to shadow invariant manifolds of more general objects (see for instance [Mar96 , 
ICre971 IFMflDl ICreOll IGZ041IZG041IB 0 IO 6 IIDGR131 ISabl.Sj ). Nevertheless, such results deal with 
invariant manifolds of objects which are (normally or partially) hyperbolic and not parabolic as 
in the present setting. 
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We generalize such result to the parabolic setting. Nevertheless, we are dealing with the 
simplest case, that is, fixed points for the 4-dimensional Poincare map, which correspond to 
periodic orbits for the flow. In the general hyperbolic setting usually one considers C 1 Lambda 
lemmas even if for the shadowing results only a C° version is needed. The reason is that to 
have C° estimates including the tangential directions, one needs to compute the C 1 estimates 
at the same time. In the present setting, since we are dealing with the invariant manifolds of 
fixed points of maps, there are not tangential directions and, therefore, one can directly compute 
the C° estimates a la Shilnikov (see Sil67j b This is a considerable simplification since in the 


parabolic setting the classical C 1 estimates are no longer true, as shown in [ GK12] , 

Lemma 3.1. Let T be a curve which transversally intersects W s { at P G W’ s (A ao y; 0 ) 
for some A QOiGq G a[ Gi ’ G2 1. Let Z be a point on W u {A a0t c 0 ). For any neighborhood U of Z in 
M 4 and any e > 0, there exists a point a € B e (P) n T and a positive integer n (which depends 
on Z, e and U) such that V n (a) G U. 

As a consequence W u (A a0t c 0 ) C Uj>oVi (T). 

The proof of Lemma 13.11 is placed in Section The lemma follows from Proposition 14.31 

Remark 3.2. By reversing time, one can consider curves T intersecting transversally W u {h( Gl ’ G2 ^) 
at P G W u (A a0t G 0 )- Then, we have the following statement. Let Z be a point on W s (A QO iGo ). 
For any neighborhood U of Z in M 4 and any e > 0 there exists a point a G B e (P ) n T and a 
positive integer n (which depends on Z, e and U ) such that V~ n (a) G U. 

As a consequence W s (A aQ! G 0 ) C U />o'P” JI (r). 

Now it only remains to apply Lemma 13.11 to shadow the transition chain given in Corol¬ 
lary [23J This argument is standard, one can see, for instance, [DdlLSOO] . We include it here 
for completeness. 

Proposition 3.3. Let {A akt c k }k>o be a family of parabolic fixed points in A[ Gl ’ G2 ] of the 
Poincare map V in (usd such that, for all k, W u (A akt c k ) intersects transversally W s (A[ Gl,G2 l) at 
Pk G W s (A afc Gj . +1 ). Consider two sequences of real numbers {5fc}fc>o an d dk,dk > 0. 

Then, there exist a G Bs 0 (A a0j G 0 ) and two sequences of natural numbers {A^j^o, {Nk}k>o, 
N k < N k < N k+ 1 < N k+ i for all k, such that V Nk (a) G B Sk (A ak:Gk ) and V^ k (a) G B^(p k ) for 
all k. 


Proof. The proof follows closely the arguments in [DdlLSOO] , 

We are going to construct a sequence of nested non-empty compact sets Ui C Bs 0 (A a0j G 0 ) 
with the following property: if p G Ui, its forward orbit by V visits the balls B$ k (A aki G k ) and 
B$ k (pk) for 0 < k < i. 

Let xo G i? 5 0 (A Q0i G 0 ) n W s (A q , 0i g 0 ). We can choose an open neighborhood Uq of xo such 
that Uq C Uq C R 5o (A QO iGo ). By Lemma 13.11 (see Remark 13. 211 , there exist a point y$ G 
B$ 0 (po) ClT s (A Ql]Gl ) and an integer no > 0 such that V^ n °(yo) G Uq. By continuity of the map 
V, there exists an open neighborhood Vo of yo such that V~ no (Vo) C Uq. Since yo G W s (A ai;Gl ), 
there exists ni > 0 such that V ni {yo) = x\ G B$ 1 (A Q1)Gl ) n W s (A Q , liGl ). By continuity, there 
exists an open neighborhood U\ of xi, such that V~ ni {Ui) C Vo- We define U\ = V~ n °~ ni {Ui). 
By construction, U\ C Uq C Uq and if x G Li, V n ° (x) G Vo C B~^(pq) and V n ° +ni (x) G U\ C 
B$ 1 (A aij G 1 )- We take Nq = 0, Nq = no and N\ = uq + n\. 

The proof follows by induction. Then, any point in C\i>oUi ^ 0, satisfies the claim. 

□ 
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Proof of Theorem 11.11 is a direct consequence of this proposition. 

Proof of Theorem 11.11 Let {A akj G k }k>i be one of the bounded orbits given in Corollary 12.51 
Since A ak ,G k = 5 + (A a i ._ 1) G fc _ 1 ), k > 1, the unstable manifold of A. ak _ lt G k _ 1 intersects transver- 
sally the stable manifold of the cylinder in a point pk in the stable fiber of A a k ,G k - 

Now we apply Proposition 13.31 We take 5k = 1/k (or any other positive sequence with limit 
0) and 5k = 5 small enough so that the balls Bg(pk) do not intersect the cylinder A^ 1 ’® 2 !. Let 
Nk and Nk be the natural numbers and a € Bg 0 (A a 0 ) c 0 ) be the point given by Proposition 13.31 
Then, the orbit of a is oscillatory. Indeed, liminffc _ ) . +00 dist('P 7 Vfe (a), A. akt Q k ) = 0 since 5k -* 0 
as k —> +oo and hmsup^+^dist^^^a), A afe )( 3 fc ) > 0 since the balls Bj(pk) do not intersect 
{x = y = 0}. Pulling back the McGehee change of variables r = 2/x 2 and considering the 
original coordinates, we obtain oscillatory orbits. □ 

4 AC 0 Lambda lemma: proof of Lemma 13.1 

4.1 Local behavior close to infinity 

System (ED can be written as r = 2 /x 2 , 

1 3 

x = ~jx y 

y = i G 2 x 6 + d r U (2aC 2 , a, t) = — + x e O± 

8 4 

a = -Gx 4 
4 

G = d a U (2x“ 2 , a, t ) = (3(a, t)x e + x s O\, 

where /3 is a function which is 2 - 7 r-periodic in its variables and Ok stands for d?(||(a;,y)|| fc ). 

To straighten the lowest order terms, we make the change 

Q = \{x-y) 

P = \{ x + y) 

and to make a a central variable, we consider the new variable 6 = a + Gy. Then, we have the 
new system 

q = j{q + pf(q + (q + p) 3 Oo) 
i> = -7 {q + p) 3 {p + {q + pfo 0 ) 

« (11) 

e = ( q + pfo 0 

G=(q+p) 6 O 0 

i = 1. 
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This system is a particular case of a system of the form 


q = \(q + p) 3 (.Q + (q + p) 3 Oo) 

p = -^{q + p) 3 (p + (q + pfo 0 ) ( 12 ) 

z = ( q+pfOo 
i = 1 


where ( q,p,z ) € R x M x K, K C M n a compact set. Now Of, = 0(\\(q,p)\\ k ) are T-periodic 
functions on t and the estimate is uniform for z G K. 

Notice that for any zo £ the set 

A Zo = {q = P = 0, z = z 0 , t G 1} 

is a periodic orbit of system (11211 . In this section we study its invariant manifolds. The set 
A = U 2oe jfA zo is invariant. In the case case of system dill) , when zq = (ao,Gx>), the set A 
corresponds to the “parabolic infinity”. 

Proposition 4.1. Consider the system m The set A 2(3 possesses invariant stable and unstable 
manifolds, Wf Q and Wf 0 . More concretely, 

= {( q,p,z,t ) | {p,z) = 7 u (q,z 0 ,t), q G [0,g 0 ),i G T} 

where 


1. 7 “ is C°° with respect to q and analytic with respect to zq and t, 

2. vr p 7 “(g, zq, t) = 0(q 2 ), TT z r y u (q, zq, t) = zo + 0(q 3 ), where n p and n z are the corresponding 
projections. 

The analogous statement holds for W s , as a graph over p. 

Proof. It is an immediate consequence of Theorem IA.71 in the Appendix. Indeed, after changing 
the sign of time and introducing the new variables 


2 = 


q+p 


(Z - Zq), 


system m becomes 


q = ~(q + pf(q + {q + p) 3 o o) 
p = \{q + pf{p + {q + pfo 0) 

5 = j(q + p) 2 (q-p)z+ {q + p) 5 Oo 

i = 1 

and depends analytically on zq and t. It satisfies the hypotheses of Theorem IA.7I with N = 4 
and c = 1/4. Then, the origin of (1131) has an invariant stable manifold parameterized by 

7 (u, zo,t) = (u, 0,0,0) + 0(u 2 ). 
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From this parametrization, we have that the invariant manifold is the graph of a function over q 
satisfying ( p,z) = 0(q 2 ). The claim follows simply restoring the original variables. 

The proof for the stable manifold is analogous. □ 

Once we have the existence and regularity of the invariant manifolds, their straightening can 
be easily accomplished. 

Proposition 4.2. There exists a C°° change of variables which transforms system m into 

q = f{q,pfq{i + Oi) 

p = -f(q-,pfp(X + Oi) 

2 = f(q,p) 2 qpO i 

t = l 


defined for q,p > 0, where f(q,p) = q + p + O 2 is a C°° function in {q,p > 0} with bounded 
derivatives. 

Proof. The change is the composition of two consecutive changes of variables. The first one is 
defined as follows. By Proposition 14.11 the map 

T : (q,z 0 ,t) ^ (q,TT z 7 u (q,z 0 ,t),t) 

is a C°° diffeomorphism from {(q,zo,t) | 0 < q < 5, t € T} onto its image, for some 5 > 0. By 
item 2 of Proposition 14.11 there exists 5' > 0 such that {(q,z,t) | 0 < q < d'} C T({(q, zo,t) \ 
0 < q < <5}). Hence, 

(q,zo(q,z,t),t) = r _1 (g,z,t) = (q,z,t) + 0{q 3 ) 


is a diffeomorphism. 

The first change that we perform is given by 


Q = q, P = P-'K p 'y u (q,z 0 (q,z,t)), z = z 0 (q,z,t). 

It is clearly a change of variables if q is small, and straightens the unstable leaves. Analogously, 
one straightens the stable leaves. The second change is the identity on the unstable manifold. □ 


4.2 The Lambda lemma 

To state the Lambda Lemma 13.1 1 it is more convenient to work with the Poincare map 


V : {t = t 0 } —» {t = t 0 + 27 t } ( 15 ) 

associated to the flow of the system (114[) . 

Lemma 13. II follows from the next statement. 

Proposition 4.3. Consider a point Q = (qf,0,zo) G W u (0, 0, zq) and a C 1 curve C parameter¬ 
ized by a(5 ) = (5,ct p (5),a z (5)), 5 G [0, <5o] for some 5 q > 0, which is transversal to {q = 0} in 
a point P = a(0) = (0,po,zo). For any /)< 1, there exists a sequence {d^}fc>i \ 0 as k —>• 00 
with 0 < 5 k < do, and an increasing sequence of {n k } k > 1 C N, n k —> +00 as k —>• +00, such that 


V nk (a(6 k )) € B p (Q). 
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Proof. We first study the associated equation m and later we show how to deduce the state¬ 
ment for the Poincare map. 

As {q = 0} and {p = 0} are invariant manifolds for system (fT4l) . the region {q > 0,p > 0} 
is invariant. In this region, by Proposition 14.21 / is strictly positive for (q,p) small enough and 
satisfies \f(q,p)\ > -K’||(?,p)||, for some K > 0 depending only on the domain. 

We rescale time by setting 4§ = /(<?>p) 3 - In the new time s , in the region {q > 0,p > 0}, 
the system becomes 


q = q{l + O l ) 
p = -p( 1 + Oi) 
z = qpO 0 


(16) 


Call d ,s the flow associated with this equation. 

Fix £o > 0 and Co > 4 such that the equation (fT6l) is well defined for 

( q,p,z,t ) eO = (0, £q} 2 x B^ 0 (z 0 ) x T C M 2 x M n x T. 


There exists K > 0 such that the terms Oo and 0\ appearing in (fT6l) satisfy \Oo\ < K and 
|Oi|<K||(g,p)|| < C 1. Choose s G (0,eo) such that e = Ke G (0,1/10). Then, for any point in 
(0,£] 2 x [Gi,G 2 ]xT 2 , 

(1 -i)q <q< (l+i)q 

-(1 +i)p <p< -(1 — e)p (17) 

—Kqp < ii < Kqp , i = 1,..., n 

where z = (zi,..., z n ). 

The points P and Q introduced in the statement of the proposition can be chosen such that 
0 < qf,Po < e/10. Fix do > 0 with do < minjdo, e/10}. Fix p > 0 and define U = B p (Q) C M 4 . 

For any 0 < 5 < 5q, there exists T* = T*(e,5) > 0 such for any time in [0,T*], the orbit of the 
point a(5) under the flow T s of system (1161) does not leave the domain D. 

Applying Gronwall estimates to the first two equations of (ED with initial condition (q(0),p(0), z(0)) 
a(5) and s G [0, T*], one obtains the following inequalities 


q , (0)e( 1 e ' )s < q(s) < < 7 ( 0 )e( 1+e ^ 

p(0)e _ ( 1+£ ^ <p(s)< p(0)e^ ( ' 1 ~ £ ' )s . 


Moreover, 

< 5Ke s/5 q(0)p(0). (19) 

It is clear that the variable that leaves first the domain D is q and therefore T* satisfies 


|z(s) - 2 ( 0 )|| < 


qpO 0 


— In ^ <T* < 

+ e 5 


1 l e 

1-^ ln I- 

1 — £ 0 


Since qp < e/10, by continuity, for 6 < do there exists S = S(5 ) < T* such that T s '(a((5),0) G 
{q = qf}. Moreover, by the Gronwall estimates in (1181) . 


— In -y) < S < 
+ £ 5 - - 




1-^ ln X 

1 — £ 0 


(20) 
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For such S', using the second inequality in (fl 8 |) we get that 


1 — £ 
1+e 


p(S) < p(0)e 1+e 1o§ q f <p( 0) (J^j 

Note that, since e <C 1, the right hand side converges to zero with 6. Now we use CHD, to bound 
|| z(S) — z( 0 )||. We have that 


II z(S) - z 0 \\ < ||z(5) - 2 ( 0 )|| + || 2 ( 0 ) - Zq || 

< 5Ke s ^ 5 q(0)p(0) + ||ck|| c i<5 

< 5Kp(0)qJ {1 - e) 5 1 5 (i-«) + ||a|| c i<5 

where we have used that q( 0 ) = 5. 

We want the final point in {q = qf} to belong to B p (Q). We choose d such that this is true. 
We need 

\p{S) ~Po\<£, and \\z(S) - zoll < j- 
Then, we need to choose 5 such that 


1 —£ 
1+e 


C° ( — ) < j and 5K\\~f\\ c oq 5 f (1 e) 5^ + 


r-l- T 


C 10 s 


<p. 


All is left to prove now is that, in the original time variable (before rescaling), that is for the 
flow associated to the equation (| 14[) . the time T needed for the transition from 7 (5) to B p (Q) 
can be chosen as a multiple of 2ir. This is achieved by taking 5 small enough and allows us to 
say that what we obtained makes sense for the Poincare map. The argument goes as follows. 
Fix the starting section p = p( 0) and the final final section q = qf. The transition time S in 
the rescaled time depends continuously on 5 and, by (1201) . satisfies that S(5) —> 00 as 5 —> 0 . 
Since the orbit never leaves the domain D, we have that dt/ds > 1/e 3 . Therefore, the transition 
time T(5) —> 00 as 5 —> 0. By continuity, it must pass through a multiple of 27r. In fact, it 
must pass through a multiple of 27r infinitely many times T(5k) with 6^ —>• 0 as stated in the 
proposition. □ 
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A Invariant manifolds of parabolic points and dependence with 
respect to parameters 

Here we present a version of Theorem 2.1 in |BFdlLM07| on the invariant manifolds of parabolic 
fixed points where we also include their dependence with respect to parameters. We need this 
statement for Proposition 14.11 which provides the proper set of coordinates in which we derive 
de Lambda Lemma. Here we deal with the analytic case, which is the one relevant in the 
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present work, unlike the setting of [BFdlLM07| . where C k maps were considered. The proof of 
the theorem we present here follows the same lines of the one in !BFdlLM07| , but is simpler and 
shorter. The main difference is that, to deal with analyticity, we work with complex domains. 

We first present the statement for maps and then we deduce the analogous result for time 
periodic flows. 

Given a local diffeomorphism F from a neighborhood of (0, 0) G R x R n to R x R n such that 
F{ 0,0) = (0,0), an open convex set V such that (0,0) G dV and r > 0, we define the local stable 
set of (0,0) with respect to V as 

Wfy = {(s, y) G V n B r I F k (x , y) G V n B r , for all k > 0, lirn F k (x , y) = (0,0)}, 

V ’ r k —^CXD 


where B r denotes the ball of radius rinRx R”. 


Theorem A.l. Let U C M x R n be a neighborhood of the origin (x,y) = (0,0), where x G R 
and y G R™, A C R p an open set and F = (F 1 , F 2 ) : U x A —> R x R n a real analytic map such 
that its derivatives up to order N > 2 at the origin are independent of the parameters A G A, 
F( 0,0, A) = 0, DF{ 0,0, A) = Id 


and 


and 


D J F(0, 0, A) = 0, 
d N F ] 


dx N 


(0,0, A) =: —c < 0, 


for 2 < j < N — 1 
d N F 2 


dx N 


(0,0, A) =0 


d N F 2 

spec ^ M _ 1 0 (0,0, A) C {z G C | Re 2 > 0}. 


dx N ~ 1 dy 

Then there exists to > 0 and a C°° map K : [0, to) x A C R x 
K(t) = (t, 0) + 0(t 2 ), analytic in (0, to) x A, and a polynomial R(t) = t — ct N + c(X)t 
c analytic in A, such that 

FoK = Ko R. 


x R n of the form 
N 1 ^+2^-1 with 


( 21 ) 


Furthermore, there exists an open convex set V, with (0,0) G dV, containing the line {x > 
0, y = 0}, r > 0 and to such that the range of K is 1F~W 


Proof. Along the proof we skip the dependence of the functions on A. 

First of all, we remark that the hypotheses of Theorem IA.1I imply those of Theorem 3.1 
in [BF04j . Hence, there exists an open convex set V, with (0,0) G dV, containing the line 
{cc >0,2/ = 0} and r > 0 such that the local stable set W~ ^ is the graph of a Lipschitz function 

if : [0, so) —> R n . The fact that the set V can be chosen independently of A follows from the 
estimates of Theorem 3.1 in BF04j . where it is proven that it contains a cone whose size only 
depends on the derivatives of F up to order N. Then, if K is the solution of (1211) given by the 
theorem, since K(t, 0) = (f, 0) + 0(f 2 ), its range is the graph of a Lipschitz function, which must 
coincide with W~ . This proves the last statement of the theorem, assuming that the previous 
ones are true. 

Now we prove the existence and properties of such a K. 

Under the current hypotheses, we can apply Lemma 3.1 of |BFdlLM07] and we obtain that for 
any k > N there exist polynomials K- k : R —> R 1+n , of degree k, and R(t) = t — ct N + c(X)t 2N 
such that 

F o K^ k {t) - K- k o R(t) = <D(t k+N ). (22) 
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It is straightforward from the proof of Lemma 3.1 that the coefficients of K- k and R are analytic 
functions of the coefficients of the Taylor expansion of F at (0, 0) up to order k. Hence, they 
depend analytically on A. 

We claim 

Lemma A.2. Let a = 1 /(JV — 1). Let A be a neighborhood of A in C p in which the map F is 
analytic. There exists ao > 0 such that for any 0 < a < ao, there exists p' > 0 such that the set 
V = {t £ C | | argf| < a, 0 < \t\ < p} satisfies R(V) C V, for any 0 < p < p' and A G A. 

Proof of Lemma \A.A We remark that, for t G V, 

. jV —1 ~2N-2\ •, (1 — ct N ~ l + ct 2N ~ 2 ) 

arg(l — ct + ct ) = -tI°g | (1 _ ctJf -i +3M _ 2 j| 

= -clt^- 1 sin ({N - 1) arg t) + 0{\t\ 2N ~ 2 ) 

= —(N - l^r- 1 arg t + \t\ N ~ 1 0(a 2 ) + 0(\t\ 2N ~ 2 ). 

If p is such that 1 — (N — 1 )cp N ~ 1 > 0, for some M > 0 we have 

| arg R(t)\ = | arg(f(l — ct N_1 + ct 2N ~ 2 ))\ 

= | argt + arg(l — cf^ -1 + c^ 2jV 2 ) | 

= |(1 - (IV - l)c|t| 7V - 1 ) argf + |t| Ar ~ 1 0(a 2 ) + 0{\t\ 2N ~ 2 )\ 

< (1 - (N - l)c|t| iV “ 1 )| argt| + + M\t\ 2N ~ 2 . 


From the last inequality we obtain that, if p satisfies, 


p N ~ 1 <{N - 1 )ca 


Ma \ 


then | argii(t)| < a, which proves the claim. □ 

We choose the constants a and p accordingly to Lemma IA.21 Taking a and p smaller if 
necessary, it is clear that there exists 0 < b < c < d such that for all t G V and any Aei, 


Rd(\t\) ■= 1*1 - d\t\ N < \R(t)\ < |t| - b\t\ N =: R b (\t\). (23) 


Notice that b and d can be chosen arbitrarily close to c taking a and p small enough. 

Next lemma describes the contraction provided by the nonlinear terms. 

Lemma A.3. Let bo, do and s > 0 such that by~ l = ab^ 1 , d^ _1 = ad -1 and bos~ a = p. Then 
there exist two sequences, (6j)j>o and (dj)j>o such that for any 0 < f3 < 1, 


for i > 0, 

and 


k = bo(l + 0(r f> )), di = do{l + 0(i~ p )), i> 0 

d i ■ i bi 

(s + i + 1)“ (s + i)"’ 

jj f bj A _ bj+i f dj A _ dj .|-i 

6 V(s + *) a y “ (* + * + !)«’ V(s + *) a / “ (s+ * + i) a ‘ 


(24) 

(25) 

(26) 
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Furthermore, the sets 


Vi = <t eV 


di+l 


< 1*1 < 


(s + i + l)“ (s + *) c 

satisfy V = Uj>oVi and R(Vi) C V l+ \. Consequently, if t G Vi, for any j > 0 


(27) 


Mi + o<r*)) < < Ml + or;» 

(s + i + j + 1)" - 1 ~ (s + i + j) a 


(28) 


Proof of Lem,ma \A.fA The two relations in (12611 define the numbers bk and dk for k > 1. To 
prove (125]) we proceed by induction. For the first step, since b < d, we have that do < bo- Then, 
using (1231) . the fact that Rd(r), Rb(r) < r, for 0 < r < ro and are strictly increasing, 


d i 

(s + l) a 



Since 

* (V) - r A_ 

\s a J (s + l)“’ 

we have that d\ < b\ and we can perform the induction procedure. 

Relations (fMl) follow from Lemma 4.4 in jBFdlLM07| . 

From pil) and (f251) follows that V = ^k>oVf ; - Inclusion R( 14) C I4+i follows from OUD and 
inequality (l24l) . □ 


We choose k > 2N — 1. Let K- = K- k : R R 1+n , polynomial of degree k given in (l22l) , 
and R(t) = t — ct N + ^(A)^ 2 ^" 1 be such that 


E{t) := F o K-(t) - K- o R(t) = 0{t k+N ). 


(29) 


Since from now on k is fixed, we skip the dependence on k of K— and E. It is worth to remark 
that if one chooses k! > k and finds the corresponding K- of degree k', its terms up to degree k 
coincide with the former (since we keep R fixed in the construction). 

We want to find a solution of 


F o (K^ + <p) — (K- + tp) o R = 0. 


We rewrite this equation as 
where 

and 


£(<r) = 

£(<p) = (DF o K-)tp — Lpo R 


?{#) = —E — F o (K^ + ip) + F o K- + (DF o K-)tp. 

The operator dehned by (Kill) is linear. A formal inverse is given by the formula 

0(VO = [ II(^)- 1 0 K~ o R*] if o RE 

j >0 \i=0 J 

(see |RFdlLMn7| l. 


(30) 

(31) 

(32) 


( 33 ) 
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As explained at the beginning of Section 4.4 in }BFdlLM07j . we make the rescaling y = by 
for some parameter 5 and from now on we work in this rescaled variable without changing the 
notation for K— and F. For any a > 0, choosing appropriately the norm in R x R n and taking 
a , p and 5 small enough, we can assume that, for some C > 0 , 

||Di ? ~ 1 (AT-(t))|| < 1 + (N + (r)d|f |' /v_1 + C\t\ N , for all f e V. (34) 


Lemma A.4. Let {V)},>o be the family of sets defined in (1271) and s > 0 the number defined in 
Lemma \A .31 There exists C > 0 such that for any l > 0, t G V), the following inequalities hold 


PJ(DF ) -1 o K- o Bfft) 


i =o 


< C 


s + l+j 
s -\-1 


(N+cr)adb 


Proof of Lemma A.f. By Lemma I A. 31 R l (Vi ) C Vj+j. By (134|) . and using the definition and 
properties of the sets V) in Lemma IA.31 if t E V), 


| (DF)- 1 o K- o i?(t)|| < 1 + (N + a)d 


(s + i + iy 


N-l 


+ c 


(s + i + iy 


N 


< 1 + (.N + a)adb 1 + 


C 


s + i + l (s + i + /) 1+7 ’ 
where 7 = min{a,/3} (see (|24l) l. Then, redefining C , 


il(DF)-\K<- mm 


i =0 


< exp I XI log ( 1 + <iV + + 


\i =0 
' 3 


c 


s + i + l (s + i + iy+7 


(.N + a)adb 


-1 


+ 


C 


s + i + l (s + i + 0 1+7 

s + j + l^(N+a)adb^ c 


1 


7 \(s+ 1 — 1)7 ( S + l + j)7 


< exp 

\j=0 

< exp log 

s + l — 1 

which implies the claim with a suitable C. 

In order to solve equation (l30l) . we introduce the space 

X m = {ip : V x A —> C 1+n | ip analytic, sup \t~ m ip(t, A)| < 00 }, 

teV,\eA 


1 


□ 


which is a Banach space with the norm 


IMI m= sup \t m (p(t,X)\. 
t£V,\eA 

By ([29]) . the function E in (|29l) belongs to Xk+N- 

Lemma A.5. Assume m G N satisfies m > N + (N + a)adb ~ 1 + 2. Let F and Q be the 
operators defined in (ED and (1331) . Then Q : A m _)_jv-i X m is linear, bounded and C o Q = Id 
On Xm+M—l- 
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Proof of Lemma \A.fft. Let ip G X m+ N-\. For any l > 0, t € V), by Lemma [A.41 and inequali¬ 
ties (1281) . we have that for some constant C independent of ip and l, 


mm i<E 

j> 0 


'J(DF)- 1 o o R\t) 


i =0 


| ip o W (t) | 


i>o 

< 


s + l + j\ {N+a)adb 1 /6 0 (l + O(^^))\ m+iV ' 1 

5 + Z ) \ (s + l+j) a ) ^Wrn+N-l 


c 


(g _|_ ^a(m+iV-l)-l IIV’llm+AT-l 
Hence, since V = U/>oV), 


(s + i + iy 


— SUp sup 1 1 m G{lp){t)\ <C sup ^ na (m+N-l)-l IIV’llm+TV-l < C\\lp\\ m+N -i ■ 
l >0 teVi,x&A l >° I s + v 


The above calculations show that the sums that define Q are absolutely convergent and can be 
reordered if ip £ X m+ N-i- Then a simple computation shows that C o Q = Id □ 

To complete the proof it is enough to check that the operator T defined in (1321) is well defined 
and Lipschitz between appropriate Banach spaces with Lipschitz constant small enough. Recall 
that J-^O) = —E € Xk+N- 

Lemma A. 6. Let r > 0 be such that the map F is analytic and bounded at B r ( 0,0) x A C 
R x R n x R p , where B r ( 0,0) is the ball of radius r centered at the origin mix R n . Let 
B C X/, + i the ball of radius R > 0. Assume that p < (r / R) l ^ k+1 ' > . Then F : B —>• X^+n is well 
defined and Lipschitz with 

lip T < CRp k ~ N+1 . 

for some C > 0 independent of R and p. 

Proof of Lemma \A. (A Let ip € B. It satisfies |<^(f)| < ||^||fc + i|t| fc+1 < Rp k+l < r. Hence, the 
function F(p) is well dehned and 

HJX^Ilfc+iV < ||-E||fc+iV + sup C\t\- k ~ N \p(t)\ 2 < Pllfc+AT + sup C'|f| /c_JV+1 ||< / ?||f +1 < oo. 

teV.XeA t£V,X&A 


As for the Lipschitz constant of J 7 , for any (p, p' € B, since k > 2 N — 1, 


- F{p')\\ k+N < sup C\t\ k N max{\p(t)\,\p’(t)\}\(p - p')(t)\ 

££ V, 

< sup C\t\ k ~ N+1 R\\ip - p'\\k+i 
tev,X£A 

<CRp k ~ N+1 \\p-p'\\ k+1 . 


We consider the fixed point equation 


p = Q o F{p). 


□ 


(35) 
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Let k > min{2A — 1, N + (N + a)adb~ 1 + 1} and choose m = k + 1 in Lemma IA.51 In 
Lemma [Oil take R = 2\\Q\\ ||£'||fc + jv and p < min{(r/i2) 1 /( fc+1 ), (R||£/||) _1 /( fe_iV+1 )}. Then, the 
map Q o F : B —> B is well defined and Lipschitz with lip (Q o J) < 1. Hence, equation (1351 has 
a unique solution p* £ B C Xk+i- By Lemma [A.51 p* is a solution of equation (1HUD . 

Up to this point we have found a solution of equation (121 |l of the form K— + ip*. Since K- 
is a polynomial of degree k, it is C k at the origin. Also, since ip* £ Xk+h lim t _>.o D^p*{t) = 0, 
for 0 < j < k. Hence, p* is also C k at the origin. Since k is arbitrary and a solution of the 
equation (1351) for k! > k also provides (conveniently rewritten) a solution of the equation for k, 
taking p smaller, if necessary, we have that K is C°°, for real t, at the origin. 

This completes the proof of Theorem IA.1I 

□ 


In what follows, given T > 0, T stands for the torus R/T. 

Given a local T-periodic vector field A from a neighborhood of {(0, 0, t), (£T}cKxK n xT 
to R x R n such that X (0,0, t) = (0,0) for all t £ T, an open convex set V C R x R n such that 
(0,0) £ dV and r > 0, we define de local stable set of (0, 0) with respect to V as 

W~ = {( x,y ) £ V n B r | (j) t T (x,y) £ V n B r , for all t > 0, lim (f t T (x,y) = (0,0), r £ T}, 

v u t—> oo 


where B r denotes the ball of radius r in R x R n and (j>t, T is the flow of the vector field X, that 
is, 

—(ft,T(x,y) = X(<t> ttT (x,y),t), <t> T , T {x,y) = {x,y). 

Theorem A.7. Let U C R x R ra be a neighborhood of the origin (x, y) = (0,0), where x £ R and 
y £ R n ; A C R p an open set and X = (A 1 , A 2 ) : UxTxd->lx R n a real analytic T -periodic 
vector field, such that its derivatives up to order N > 2 are independent of the parameters A € A 
and t at (x, y) = (0,0), A(0,0, t, A) = 0 


D J A(0,0, t, A) = 0, for 1 < j < N — 1 


and 


and 


d N X k 

dx N 


(0,0,t,A) 


—c < 0, 


d N X 2 

dx N 


(0,0, t, A) = 0, 


spec 


qN 

ftfffFTXofj (O’ 0, t, A) C {z € C | Rez > 0}. 


Then there exists an open convex set V, with (0,0) £ dV, containing the line {x > 0,y = 0} ; 
r > 0 and a C°° map K : [0, so) xTxiclxIxP-xRx R n of the form K(s, r, A) = 
(s,0) +0(s 2 ), analytic in (0, sq) x A such that the range of K is 


Proof. Let F T be the Poincare map of the vector field A associated to the section {t = t}. It 
depends analytically on r and A. A simple computation shows that it satisfies the hypotheses 
of Theorem IA.11 Let K(s,X,t ) be the solution of the invariance equation 


F t o K(s , A, r) = K(R(s , A, t), A, r) 


provided by Theorem lA.il Its range is precisely W~ 


□ 
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B Formulas for the scattering maps of the circular problem 


In this section we obtain formulas for the scattering map of the circular problem. As the 
scattering maps 5^ are defined in the extended phase space and send s to s iI)dlLS08| , they can 
be restricted to the section £ (see (jS}) giving rise to symplectic maps on the cylinder a[ G i,+0 °]. 
Recall that the result in [Ddll.SOSl deals with a normally hyperbolic cylinder, whereas in our case 
the manifold A is a “normally parabolic invariant cylinder” with stable and unstable manifolds. 
The arguments in |DKdlR,S14] extend the proof to this case. Moreover, as it is explained in 
jDKdlRS14j . the conservation of the Jacobi constant implies that G is a first integral of the two 
scattering maps. 

Since the scattering maps are symplectic and G is preserved, they must be of the form (|9|). 
Thus, it only remains to obtain formulas for the functions / . To compute these functions, we 
use [DdlLSOB] , Note that this paper provides formulas for the scattering map in terms of the 
Poincare potential in a regular perturbation regime. The regular perturbative approach used 
in |DKdlR S14| is to consider system fl2]) as an order y perturbation of the two body problem. As 
has been explained in Section [TJ here we use a different nearly integrable regime, also considered 
in [QMS 15] for the circular problem, which allow us to deal with arbitrary values of y. We 
proceed as follows. 

Recall that we already know that the scattering map is defined for a € T and G > G\. 
To provide formulas for f± we need to rescale the variables to transform the problem into a 
perturbative setting. Consider the scaling 

x = Gj" 1 ®, y = G'j~ 1 y, a = a, G = GiG 

and the change of time s = Gft. It transforms system ([2]) for eo = 0 into a system of the same 
form with a new potential 

x 2 ~ 

U(x, a, t\Gi) = G\U(G^[ l x, a, Gft; 0) = — + A U(x, a, t\ G±) 


where A U(x, a, t;G i) = G\AU(G l l x , a, i; G\) and is given by 

/ 

1 — it 

+ 


x 


AG(x,a,t;Gi) = — 




V (l - G\ x2 cos ^ + ( 2 % x2 ) 2 ) ^ (l + Tj#® 2 cos 0 + X 2 ) 2 ) 


1/2 


- 1 


with (j> = a — G\t. Expanding A U in G^ 1 one can see that it can be written as A U = G^ 4 AV 
where AV is bounded uniformly on Gi —>■ + 00 . Therefore, system ([2]) with eo = 0 in these 
rescaled variables can be seen as a G(Gj” 4 )-perturbation of the two body problem. Nevertheless, 
this limit is singular since AG, and therefore AV, is 2ir/G \-periodic in time. When Gi —> + 00 , 
the frequency of the perturbation blows up. Hence, we need to adapt the theory developed in 
[DdlLSOS] to our setting, along the lines in |DdlLS00 l 

To transform this problem into a regular perturbation one, we introduce 6 = Gj~ 4 , which we 
consider as an independent parameter. Later we will recover the true value of 6. We write the 
potential U as 

x 2 

y+ iAV 
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and, then, system (J21) is a 0(<5)-perturbation, 27r/Gf-periodic in time, of an integrable sys¬ 
tem with a “normally parabolic” invariant manifold Ai Gl,+0 °) with a 3-dimensional homoclinic 
manifold that can be parameterized by the homoclinic orbits 

(xh(i; Go),yh{t\ Go), ah(i; ao, Go), Gh(i; Go), sq + Gft), 


for any Go, ao, so- Here (xh, yh, «h, Gh, Sh) is the homoclinic orbit given in Lemma [2721 and we 
have changed the last component to take into account the period of the perturbation. 

Since now system (J2J) is a 5-regular perturbation of an integrable system with a homoclinic 
manifold, we can use the perturbative arguments (using deformation theory) in Theorem 32 
of fDdlLSOS] . We consider the Poincare function associated to the homoclinic manifold 

L (a 0 , Go, sq, a + cr; Go), «h(t + cr; ao, Go), so + Gftj dt 

= G\ f AU (x b (t + a;Go),a h (t + a-,a 0 ,Go),s 0 + Gft\ dt 

J — OO ' 

— G\ f A IA (g± ^Xh(t + cr; GiGo), Go), so + dt 

= G iJ AU ^Gj'^h ( r "g 3 1(J ; ^ 0 ) ,ah ( r ~ I ^ 1<T ’ 5o ’^°) ’^0 + r ) dr 
= G\ J Au(x h (r + Gla;GiGo^,a h (r + Gla-,a 0 ,GiGo^,so + r^ dr 

= G 3 L ^ao, GiGq, sq, G 3 cr; 0^ = G\L ^ao, G'iGq, sq — G 3 cr, 0; 0^ , 


where we have used that the formulas in Lemma 12.21 imply 



and L is the function introduced in ([5|). The corresponding reduced Poincare functions are 
obtained by evaluating the function a L( ao, G,s, o , cr) at its nondegenerate critical points 0 ±, 
given by 'sq — oiq — G\a*_ =0 and sq — Sq + 7r — G\= 0 (see (]6|1). They are defined as 


C*_(ao, Go) — L(ao, Go, so, G 1 3 (so — ao)) — G\L{a o, GiGo, «o, 0; 0) 

C* + (ao, Go) = L(a o, Go, sq, G l 3 (7t + sq — ao)) = G' 3 L(ao, GiGq, ao — vr, 0; 0). 


Since L, and, consequently, L, only depend on the angles through ao — so + cr, C± only depend 
on GiGo- 

r (GiG 0 ) = G 3 L(0,GiG 0 ,0,0;0) 

G;(GiGo) = G 3 L(vr,GiGo,0,0;0). 

Then, the generating function of the scattering map for the rescaled system is given by 

St sc (a 0 ,Go) = Go 5 0 + <5*4 (Gi G 0 ) + 0(5 2 ). (36) 

Using the results in |GMSl5j . we have that 

C*±(GiGo) = G 3 L 0 (GiG 0 ) ± G 3 Li,_i(G!Go) + O (g^Go)-! e - 2i£l ^) 
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where 


Lq(GiGq) — n( 1 — fi)- 


ir 


2(GiG 0 ) 3 

£i,-i(GiGo) = /r(l — /i)( 1 — 2p) 


l + O{(GiG 0 ) 


-4 


7T (GjGfrjf 


l + O((G!G 0 ) 


-l 


f 8(GiG 0 ) 

Using (1361) . we have that the rescaled scattering map of the circular problem is of the form 

c± (& o' 


Go 


a o + /^(G o) 

, Go , 


with 

7 ± (Go) = %/i(G!Go) + O (<5 2 ) = -Sfi( 1 - + o (,5Gr 4 G- 8 ) + G (<5 2 ) . 

where O refers to uniform bounds Go £ [1, +oo). 

We recover the value of <5 = Gf 4 and we scale back the variables. Then, we obtain formula Q 
with 

Q fTT 

/ ± (g) = -m(i - /i)^4 + ° ( G r 8 ) • 

To recover the claim of Proposition I2.3L we need to replace the error term O (G^f 8 ) by O (G~ s ). 
To do so, we consider the above procedure restricted to compact sets. That is, for any k € N, 
we choose G £ [kG*, (k + 1)G*] and thus we take 5k = G^ A = ( kG*)~ 4 . Then, for any G £ 
[kG*, (k + 1)G*], 5k < CG ~ 4 , for some G > 0 independent of k. 
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